home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DSIDI.z / DSIDI
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSIIIIDDDDIIII((((3333FFFF))))                                                            DDDDSSSSIIIIDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSIDI   - DSIDI computes the determinant, inertia and inverse of a double
  10.      precision symmetric matrix using the factors from DSIFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE DSIDI(A,LDA,N,KPVT,DET,INERT,WORK,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAA DOUBLE PRECISION(LDA,N)
  20.         the output from DSIFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array A.
  24.  
  25.      NNNN INTEGER
  26.         the order of the matrix A.
  27.  
  28.      KKKKPPPPVVVVTTTT INTEGER(N)
  29.         the pivot vector from DSIFA.
  30.  
  31.      WWWWOOOORRRRKKKK DOUBLE PRECISION(N)
  32.         work vector.  Contents destroyed.
  33.  
  34.      JJJJOOOOBBBB INTEGER
  35.         JOB has the decimal expansion  ABC  where
  36.         if  C .NE. 0, the inverse is computed,
  37.         if  B .NE. 0, the determinant is computed,
  38.         if  A .NE. 0, the inertia is computed.
  39.         For example, JOB = 111  gives all three.  On Return Variables not
  40.      requested by JOB are not used.
  41.  
  42.      AAAA contains the upper triangle of the inverse of
  43.         the original matrix.  The strict lower triangle
  44.         is never referenced.
  45.  
  46.      DDDDEEEETTTT DOUBLE PRECISION(2)
  47.         determinant of original matrix.
  48.         DETERMINANT = DET(1) * 10.0**DET(2)
  49.         with 1.0 .LE. DABS(DET(1)) .LT. 10.0
  50.         or DET(1) = 0.0.
  51.  
  52.      IIIINNNNEEEERRRRTTTT INTEGER(3)
  53.         the inertia of the original matrix.
  54.         INERT(1)  =  number of positive eigenvalues.
  55.         INERT(2)  =  number of negative eigenvalues.
  56.         INERT(3)  =  number of zero eigenvalues.  Error Condition
  57.  
  58.      AAAA division by zero may occur if the inverse is requested and  DSICO  has
  59.      set RCOND .EQ. 0.0 or  DSIFA  has set  INFO .NE. 0 .  LINPACK.  This
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSIIIIDDDDIIII((((3333FFFF))))                                                            DDDDSSSSIIIIDDDDIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  75.      Nat. Lab Subroutines and Functions BLAS DAXPY,DCOPY,DDOT,DSWAP Fortran
  76.      DABS,IABS,MOD
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.